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Abstract 

o : 

The whole Hilbert state space of an n— qubit spin system can be divided 
^ . into n + 1 state subspaces according to the angular momentum theory of 

quantum mechanics. Here it is shown that any unknown state in such a 
^ | state subspace, whose dimensional size is proportional to either a polynomial 

; or exponential function of the qubit number n, can be transferred efficiently 

into a larger subspace with a dimensional size generally proportional to an 
q ■ exponential function of the qubit number by the multiple-quantum unitary 

^-i ! transformation with a subspace-selective multiple-quantum unitary operator. 

The efficient quantum circuits for the subspace-selective multiple-quantum 
O \ unitary operators are really constructed. 

^ : 

1. Multiple-quantum transition between state subspaces of the 
Hilbert space 

cd Multiple-quantum transition processes [1] are closely related to the quan- 

Q-f turn computation [2, 3]. In a quantum system such as a coupled multiple- 

particle two-level quantum system, e.g., an n— qubit spin system which may 



consist of n coupled nuclei of spin 1/2, quantum state transfer or transition 
between different states of the Hilbert space of the system generally consists 
of a variety of multiple-quantum transition processes except those nonco- 
herent transition processes which are usually non-unitary processes [1]. A 
quantum search algorithm [4] running in a quantum system usually starts 
at an initial state such as a superposition, then performs a sequence of 
known unitary operations and the oracle unitary operation to convert the 
initial state into the marked state in a high efficiency, and finally makes a 
quantum-mechanically measure to output the marked state. A variety of 
quantum state transition or transfer from one state to another really occur 
in the quantum system during the quantum search process. The quantum 
search problem usually takes the whole Hilbert space of a quantum system 
as the search state space in which the unknown marked state is searched 
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for. The computational complexity of the quantum search problem is closely 
related to the complexity of the multiple-quantum transition processes be- 
tween different subspaces of the Hilbert space of the n— qubit spin system. 
The quantum state transition or transfer generally could not be efficient for 
an unknown state, e.g., the marked state in the quantum search problem 
from a state subspace of the Hilbert space whose dimensional size is pro- 
portional to an exponential function of qubit number n of the system into 
a smaller state subspace whose dimensional size increases polynomially as 
the qubit number, because a quantum state in a known subspace with a 
polynomial dimensional size can be determined in polynomial time. The 
multiple-quantum unitary transformation has been used to help design of 
quantum circuit [5] and quantum algorithms in quantum computation [2, 3]. 
The multiple-quantum spectra could be used to output results of quantum 
computation [3]. Using multiple-quantum unitary transformation the quan- 
tum search space of the quantum search problem [4] could be reduced from 
the whole Hilbert space of the n— qubit spin system to its largest subspace 
whose dimensional size is still much smaller than that of the whole Hilbert 
space, and hence this reduction could speed up the quantum search process. 
The more important is that the multiple-quantum unitary transformation 
provides a useful method to manipulate quantum state transition or transfer 
from a state to another or from a subspace to another in the Hilbert space 
in quantum computation. 

The Hilbert space of an n— qubit spin system has a dimensional size 2 n 
which increases exponentially as the qubit number n. It has 2 n conventional 
computational bases that can be used to represent 2™ numbers or elements 
in quantum computation. This is a large search space in the quantum search 
problem. The quantum search space will be reduced if the whole Hilbert 
space can be divided into small subspaces. In quantum computation this 
subspace reduction for the Hilbert space has been proposed to help design 
for the quantum search algorithm [2, 6]. There are a number of symmetric 
property of a quantum system to achieve the decomposition for its Hilbert 
space. The important one is the rotation symmetry in space of a quantum 
system [7-10]. The rotation symmetry in spin space of an n— qubit spin 
system which consists of n spins— 1/2 may be used to guide the decomposi- 
tion for the Hilbert space into its small subspaces. The angular momentum 
theory of quantum mechanics gives a detailed description for the rotation 
symmetry in space of a quantum system [7-10], and according to the angular 
momentum theory the whole Hilbert space of the n— qubit spin system may 
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be divided into (n + 1) state subspaces, and each of which may be formed 
from the complete set of the eigenstates (l^^)} of the total magnetic quan- 
tum operator I z = Ylk=i ^ kz w ith a common eigenvalue M which satisfy the 
eigen-equation: 

I z \q M ) = M\q M ), (H=l). (1) 

The total magnetic quantum number (or the eigenvalue) M has (n + 1) differ- 
ent values and can take n/2,n/2 — 1, , —n/2, —n/2 for the n— qubit spin 

system. Each value of the quantum number M remarks a state subspace. 
The state subspace with the total magnetic quantum number M = n/2 — k 
is denoted as S(M = n/2 — k) or simply as S zq (k) with k — 0, 1, n. Since all 
the states in a subspace take the same value of the total magnetic quantum 
number there is not a change for the value of the total magnetic quantum 
number in quantum transition between different states within the subspace 
and hence the quantum transition is a zero-quantum transition. A quan- 
tum state transition from a subspace to another is a nonzero-order quantum 
transition since the total magnetic quantum number value is changed when 
a state is transferred from a subspace into another. According to the an- 
gular momentum theory the dimensional size for the subspace S zq {k) with 

k — 0, 1, ...,n is given by ^ ^ ^ for the n— qubit spin system, which is de- 
noted as d(M = n/2 — k) or simply as d{k). Among the (n + 1) subspaces the 
two smallest subspaces are S zq (0) and S zq (n), whose dimension is one. The 
next smallest subspaces are S zq (l) and S zq (n — 1) which have the same dimen- 
sional size n. For a spin system with an even qubit number n the largest sub- 
space is S zq (n/2) and its dimensional size is d(n/2) = n!/[(n/2)!] 2 . For a spin 
system with an odd qubit number the two largest subspaces are S zq ((n— 1)/2) 
and S zq ((n + l)/2), respectively. Both the two subspaces have the same di- 
mensional size equal d((n — l)/2) = n\/{[(n — l)/2]![(n + l)/2]!}. For a large 
n the number d(n/2) « d((n — l)/2) ~ 2™/ yJim/2 by the Starling's formula. 
Therefore, the dimensional size for the largest subspace increases exponen- 
tially as the qubit number n. Since any computational base of the Hilbert 
space of the spin system can be only in one of these (n+1) subspaces the 
quantum search space now may be limited to such a subspace in which the 
marked state is. When the marked state is in those smallest subspaces such 
as S zq (0), S zq (l), S zq (2), etc., whose dimensional size increases polynomi- 
ally as the qubit number, it may be found in polynomial time. However, 
if the marked state is in the largest subspaces then it can be found by the 
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Grover quantum search algorithm [4] but this need take an exponential time 
proportional to a/ d(n/2). Therefore, this shows indirectly that an unknown 
state usually could not be transferred efficiently from one large subspace 
whose dimensional size increases exponentially as the qubit number into a 
small subspace with a dimensional size proportional to a polynomial func- 
tion of the qubit number. Actually, this quantum-state transfer is closely 
related to the computational complexity of the quantum search problem [4], 
that is, the quantum-state transfer is as hard as the latter one. However, 
the inverse quantum-state transfer process could be efficient, that is, an un- 
known state in a small subspace whose dimensional size may be proportional 
to an exponential or polynomial function of the qubit number could be ef- 
ficiently transferred into a larger subspace by a multiple-quantum unitary 
transformation. These quantum-state transfer processes involve in quantum- 
state multiple-quantum transitions between different subspaces. It seems to 
see clearly that the quantum-state transfer could be efficient if an unknown 
state is initially in those smallest subspaces such as S zq (l), S zq (2), etc. How- 
ever, it is not yet clear whether the quantum-state transfer is efficient or not 
if an unknown state is initially in those small subspaces whose dimensional 
size is also proportional to an exponential function of the qubit number. In 
the following a subspace-selective multiple-quantum unitary operator is con- 
structed that transfers completely an unknown state from a subspace into a 
larger subspace of the Hilbert space. 

Suppose that an unknown state \ty s ) is in a state subspace S(M S ) and 
needs to be transferred to a larger subspace S(M s +p) by a multiple-quantum 
unitary transformation. The state |\& s ) can be expanded conveniently in 
terms of the usual computational basis {\(p k (M s ))} of the subspace S(M S ) : 

d(Af s )-l 

|*.)= °>sk\<P k (M s )). (2) 

k=0 

If there is a p— quantum unitary operator that can convert simultaneously 
all the computational bases {\ip k (M s ))} in the unknown state \^ s ) from the 
subspace S(M S ) into the subspace S(M S + p) then this p— quantum uni- 
tary operator also can convert the unknown state |\& a ) (2) from the sub- 
space S(M S ) into the subspace S(M S + p). It is possible to construct such 
a p— quantum unitary operator. Now dimensional sizes for the subspaces 
S(M S ) and S(M S + p) are d(M s ) and d(M s + p), respectively, and hence 
there are d(M s ) and d(M s + p) computational bases in the subspaces S(M S ) 
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and S(M S + p), respectively. By using the computational bases of the two 
subspaces one can generate a subspace-selective p— quantum unitary oper- 
ator that converts the unknown state \^ s ) from the subspace S(M S ) into 
S(M S + p). Because d(M s +p) is greater than d(M s ) one can choose properly 
any d(M s ) bases among the d(M s + p) bases of the subspace S(M S +p). By 
combining these d(M s ) computational bases from the subspace S(M S + p) 
with all the d(M s ) computational bases of the subspace S(M S ) one can build 
up d(M s ) state-selective p— quantum Hermitian operators. For every pair of 
the computational bases \(p k (M s )) of the subspace S(M S ) and \ip k (M s + p)) 
of the subspace S(M S + p) a state-selective p— quantum Hermitian operator 
is built up: 

Q psk = ^(\MM s ))(MM s +p)\ + \<p k (M,+p))(<p k (M a )\). (3) 

Then the corresponding state-selective p— quantum unitary operator is gen- 
erated from the Hermitian operator Q psk by 

U psk {9) = exp(-i6Q psk ). (4) 

This state-selective p— quantum unitary operator is only applied to both the 
two states \<p k (M s )) of the subspace S(M S ) and \<p k (M s +p)) of the subspace 
S(M S + p). There are d(M s ) such state-selective p— quantum unitary opera- 
tors. All these state-selective p— quantum unitary operators are commutable 
with each other because all the computational bases are orthogonal to each 
other. Then by multiplying all these d(M s ) state-selective p— quantum uni- 
tary operators a subspace-selective p— order quantum unitary operator is 
obtained: 

d(M s )-l d(M s )-l 

U ps {6)= n U psk {6) = exp(-iO Qp°*)- ( 5 ) 

fc=0 k=0 

This p— quantum unitary operator is selectively applied to both the two sub- 
spaces S(M S ) and S(M s +p). Since the dimensional size d(M s ) of the subspace 
S(M S ) may increase exponentially as the qubit number the unitary operator 
U ps (9) (5) may contain exponentially many state-selective p— quantum uni- 
tary operators (4). First it can be proved that the basis state \ip k (M s )) of the 
subspace S(M S ) can be converted completely into the basis state \ip k (M s +p)) 
of the subspace S(M s +p) by the state-selective p— quantum unitary operator 
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Up S k{0). By expanding the unitary operator U psk {6) (4): 

U psk (9) = E + (-l + cos±9)(\<p k (M s )){<p k (M s )\ 

+\<p k (M s +p))(<p k (M s +p)\) -i2Q psk sin -0 (6) 

and using the orthogonormal condition for the usual computational bases 
one easily obtains 

Therefore, the basis state \(p k (M s )) is transferred completely to the basis state 
\ip k (M s + p)) in additional to a phase factor —i by the p— quantum unitary 
operator U psk (ir). However, any other basis state of the subspace S(M S ) 
keeps unchanged under the p— quantum unitary operator. This result can 
be further used to prove that any computational basis state \tp s ,(M s )) of the 
subspace S(M S ) can be converted completely into the basis state \ip s ,(M s +p)) 
of the subspace S(M S + p) by the subspace-selective p— quantum unitary 
operator U ps (ir) in (5). Because all the state-selective p— quantum unitary 
operators Q psk (k = 0, 1, ...,d(M s ) — 1) are commutative to each other one 
therefore has for any basis state \(p s ,(M s )) of the subspace S(M S ), by using 
(7), 

d(Af s )-l 

U ps (ir)\cp s ,(M s )) = Yl exv(-inQ psk )\cp s ,(M s )) 

= exp(-i7rQ pss >)\ip s ,(M s )) 

= -i\<p s ,(M s +p)). (8) 

Thus, the basis state \ip s ,(M s )) is completely transferred into the basis state 
\(p 3 ,(M s + p)) of the subspace S(M S + p) in addition to a total phase fac- 
tor —i by the p— quantum unitary operator U ps (n). This also indicates that 
all the basis states of the subspace S(M S ) in the unknown state \^ s ) (2) 
can be simultaneously converted completely into the subspace S(M S + p) 
by the p— quantum unitary operator U ps (tt). Then an arbitrary state of the 
subspace S(M S ) also including |\& s ) with the expansion (2) can be also trans- 
ferred completely into the subspace S(M S + p) by the p— quantum unitary 
operator U ps {it). 
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Generally, it is hard to construct an efficient subspace-selective multiple- 
quantum unitary operator that applies selectively on both two subspaces with 
exponentially many basis states. But since the larger subspace S(M s +p) has 
a larger dimensional size than the subspace S(M S ) then it could be possible 
to construct such an efficient subspace-selective multiple-quantum unitary 
operator as U ps (9) (5) by suitably choosing the d(M s ) bases {\(p k (M s +p))} 
among all the d(M s + p) (> d(M s )) computational bases of the subspace 
S(M S + p). An important subspace-selective p— quantum unitary operator 
of (5) is related to both a small subspace S zq {k) (k 7^ n/2) and the largest 
subspace S zq (n/2) of the Hilbert space of an n— qubit spin system. If the 
marked state in the quantum search problem belongs to a state subspace 
S zq {k) (k 7^ n/2) then using the known p— quantum unitary operator U ps (ir) 
one may transfer efficiently it from the subspace S zq (k) into the largest sub- 
space S zq (n/2). This really reduces the whole Hilbert space to its largest 
subspace S zq (n/2) as the search space of the quantum search problem. Be- 
low it is shown how an efficient quantum circuit for the subspace-selective 
p— quantum unitary operator (5) can be constructed in an n— qubit spin sys- 
tem. 

2. The basic unitary operations 

The basic unitary operations U (9) used to build up the subspace-selective 
multiple-quantum unitary operators (5) may be generated from the Hermi- 
tian product operators Q by the exponential mapping: U{9) = exp(—i9Q). 
The Hermitian product operators are the tension product operators of single- 
spin operators in an n— qubit spin system and generally written as Q = 
Hi H 2 ••• H n . The Hermitian operator H k is the single-spin operator 
of the fcth spin of the spin system and can be generally expressed as a linear 
combination of the single-spin magnetization operators I kx , Ik y , and Ik z as 
well as the unity operator E k : H k = a k0 E k + a kx I kx + a ky I ky + a kz I kz , while 
the single-spin magnetization operators are related to Pauli's spin operators 
a kfi by Ikfj, = \vkn (f 1 — x iUi z )- Every single-spin operator H k can be di- 
agonalized: H k = U k H k U k , and its diagonal operator is given generally by 
H k = a k (\E k + I kz ) + j3 k (\E k — I kz ). Denote Q as the diagonal operator 
of the product operator Q. Then the diagonal operator is generally written 
as Q = Hi H 2 ••• H n . The complete operator set for the diagonal 
operators Q forms the LOMSO operator subspace [11]. The complete ba- 
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sis operator set for the LOMSO operator subspace usually may be chosen 
conveniently as the longitudinal magnetization and spin order operator set 

["]:. 

Qa = {E, hz, Zhzliz, ^hzhz Imzi •••) 2 Ilzl2z---Inz\ ■ 
In addition to the longitudinal magnetization and spin order operator set 
there are also other equivalent complete basis operator sets, of which a par- 
ticularly important basis operator set is given by [2] 

Q B = {E,Df\D^ k \... } D? lh2 - k "}, 
where the diagonal operator D klk2 - hm (m = 1, 2, n; I = 0, 1, 2 m - 1) is 
defined by 

E>? lk2 - hm = (\E kl +a l ki I klZ )^E k2 + a l k2 I k2Z ) 

®-®{\E km + a l k J kmZ ), (9) 

where the indices k±, k 2 , ■ k n and k,l,m,..., are series number of spins in 
the spin system and usually are ordered: 1 < k\ < k 2 < ... < k n < n, 
1 < k < I < ... < n\ {a l k } is a quantum-state unit-number vector, a l k = ±1; 
and all the unity operator components {E k } in the product operators are 
omitted for convenience, for example, the full expressions for the product 
operators 2I kz I lz and Df lk2 should be given respectively by 

2I kz I lz = 2E 1 ®...®E k _ 1 ®I kz ® E k+1 (g) . . . 

(g) £?,_! (g) l lz (g) E l+1 (g) ... (g) E n , 

D k ^ = E 1 <S) -<S)E kl -i ®{\E kl + a l ki I klZ )®E kl+1 (g ... 
(g ^fc 2 -i ®{\E k2 + a l k J k2Z ) (g) E k2+1 (g ... (g E n . 

m 

In particular, denote that the diagonal operator D m i = Dj 2 - m =(g) {\E k + 

k=l 

a k Ikz) { m — 1, 2, n; / = 0, 1, 2 m — 1), and there is a simple denotation 
Di = D n j used also in previous papers [2, 6]. The diagonal operator D m ^ 
is called the quantum-state diagonal operator since a conventional computa- 
tional basis state of an m— qubit spin system can be characterized completely 
by the quantum-state unit-number vector {a l k } or equivalently by the diag- 
onal operator D mt i [2]. 

It is clear that both the basis operator sets Qa and Qb are equivalent to 
each other. By expanding the diagonal operator D^ lk ' 2 --- km (9) it can be seen 
that the operator is a linear combination of the product operators of the set 
Qa- On the other hand, using the operator identity: 2I kz = (\E k + I kz ) — 
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{\E k — I kz ) every product operator in the set Qa also can be expressed as a 
linear combination of the basis operators of the set Qb- 

It is known that elementary propagators built up with the basis oper- 
ators of the operator set Qa can be implemented efficiently For the el- 
ementary propagators there is a simple recursive relation for decomposing 
a multi-body elementary propagator built up with a multi-body interac- 
tion basis operator, e.g., 2 m ~ 1 Ik lZ Ik 2Z ---h m z, into a sequence of one-, and 
two-body elementary propagators. Generally, the elementary propagator 
Rk 1 k 2 ...k m {d) = exp(-i92 m - 1 I klZ Ik 2 z---h m z) (2 < m < n) can be simply de- 
composed as [5] 

Rk 1 h a ...k m (0) = exp(-i^I km _ 1 x)exp(-iirI km _ lg I kmg ) 

7T 

x exp(-i-/ fcm _ iy ) exp( y -i92 m - 2 I klZ I k2Z ...I km - 1 z) 

x exp(i^I km _ iy ) exp(iirI km _ lz I kmZ ) exp(i|j fcm _ ix ). (10) 

This recursive relation assures that any elementary propagator built up with 
the basis operator of the operator set Qa can be efficiently decomposed into 
a sequence of one- and two-qubit quantum gates. 

The basic unitary operations built up with the basis operators of the 
operator set Qb also can be implemented efficiently. It is based on the 
fact that the selective rotation operation applying only to a given state of 
any m— qubit subsystem (1 < m < n) of an n— qubit spin system can be 
implemented efficiently. It is known that the selective rotation operation 

n 

Co(0) = exp(—i9D ) built up with the basis operator D = ® {\E k + hz) 

k=l 

of the operator set Qb can be performed efficiently in an n— qubit spin sys- 
tem [2, 4]. This is an n— qubit selective rotation operation applying only 
to the known state |0i) |0 2 ) ...|0„) of the n— qubit spin system. Generally, an 
m— qubit selective rotation operation with 1 < m < n can be also performed 
efficiently in the n— qubit spin system. This m— qubit selective rotation op- 
eration is only applied to the known state |0i)|0 2 )...|0 m ) of the m— qubit 
subsystem consisting of the first m spins of the n— qubit spin system when it 
is applied to any basis state |0i)|0 2 )...|0 m )|</? m+1 )...|</? n ) of the n— qubit spin 
system (\(fi m+1 ), |^ m+2 )> \fn) take |0) or |1)). It is defined by [2] 

C m , (9)=exp(-i9D mfl ). (11) 
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There is a simple denotation: Cj(0) = C n i(9) also used in previous papers 
[2,6]. 

The efficient implementation for the m— qubit selective rotation operation 
C mt o(0) may use conveniently the reversible AND operations and conditional 
phase shift operations. The classical AND operation is really irreversible and 
it needs to be changed to the reversible one by the Bennett's method [12] 
when it is used to construct these selective rotation operations. 

More generally, any m— qubit selective rotation operation c klk2 '" km (9) of 
any m— qubit subsystem (1 < m < n) of the n— qubit spin system also can 
be performed efficiently, 

Cf lka - km (6) =exp(-i9Df lk2 - km ), (12) 

where the diagonal operator D klk2 - km (1 < m < n, < I < 2 m - 1) is a 
basis operator of the operator set Qb given by definition (9). In particular, 
C m ,i(0) = Cl 2 - m (9). The m-qubit selective rotation operation C klk2 - km (9) 
is only applied to any state |0fc 1 )...|0fc 2 )...|0fc m ) with the known state |0) of 
the kith, k 2 th., k m th spins of the n— qubit spin system. The diagonal 
operator ]j k ^- k ™ can ^ e rea jjy converted efficiently into the diagonal oper- 
ator D m . This can be achieved with the help of the zero-quantum unitary 
transformation: 

h z = V kl (n)I lz V k M + , (13) 
where the zero-quantum unitary operator is given by 

V kl (9) = exp[-i9-(2I kx I ly - 2I ky I lx )}, 

which can be decomposed into a sequence of two two-qubit quantum gates: 

V ki (9) = exp(-i9I kx Ii y ) exp(i9I ky Ii x ). (14) 

Assume that the indices in the diagonal operator jy^ k ^- km satisfy 1 < 
k\ < k 2 < ... < k m < n. If the index k\ ^ 1 then making the zero- 
quantum unitary transformation (13) with the zero-quantum unitary opera- 
tor V k i(ir) (14) with the indices k = 1 and I — k± on the diagonal operator 
D kik a ...k m the operator £)Ji*a-fc m is converted into D l k2 - km , and if k x = 1 then 
£)fcifc 2 -fcm _ jjik 2 ...k m an( j nQ un jtary transformation is needed to apply on 
the operator D^ lk2 - km . In a similar way, if the index k 2 ^ 2 then D^-^ 
is converted into D^ k ^- km by V k i(n) with the indices k = 2 and / = k 2 . 
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By making the zero-quantum unitary transformation (13) on the diagonal 
operator D klk2 --- km at most m times one can convert completely the diag- 
onal operator £> klk2 --- km into the diagonal operator .Dq 2 ---™ (= D m> o). Then 
the selective rotation operations CQ lk2 '" km (9) (12) can be implemented effi- 
ciently since C m p(6) can be performed efficiently Finally, using m pulses 
exp(— ml kx ) or exp(— iirl ky ) at most one can convert the diagonal opera- 
tor D m>l into D mfi and D klk2 - km into D klk2 - km for any index / ^ 0. This 
is based on the fact that (\E k + I kz ) = exp(~iirI kfl )(^E k - I kz ) exp(iirl kfi ) 
(fi = x,y). Then any m— qubit selective rotation operation C klk2 '" km (9) (12) 
can be performed efficiently. 

Therefore, the basic quantum gates used to construct quantum circuits 
for multiple-quantum unitary operators include three types: 

(i) the m— qubit selective rotation operations: C m fi(9) (m = 1,2, ...,n), 

(ii) the two-qubit quantum gates: exp(—iJ k i2I kz Ii z ) (k,l = l,2,...,n), 
(Hi) the single-qubit gates: exp(—i9 kjJi I kjJ ) (// = x,y,z; k = 1,2, ...,n). 
Obviously, the three types of basic quantum gates form the universal quan- 
tum gate set in quantum computation. 

Once the elementary propagators formed from the basis operators of the 
product operator sets Qa and Qb can be implemented efficiently then these 
basic unitary operations Up (9) = exp(—i9Qp) ((5 = a,b,c) built up with the 
following three types of the basic product operators also can be efficiently 
implemented: 

(i) Q a = 2'- 1 / fclMl 4 2M2 ... h m E h+1 ... E kn - 

(ii) Q b = (\E kl + alJ kllXl ) ®(\E k2 + 4 2 I k2fl2 ) ... 
®(\E kl + 4/^) E kl+1 (8) ... (8) E kn ; 

(Hi) Q c = 2 l 1 7 fclMl ... h m <g)(lE k[+1 + a{ i+i I kl+m+1 ) ... 

®(±E krn + a t k J krn t l J®E krn+1 ®...®E kn , 
where /i Q = x,y,z (a = l,2,...,n). This is because each of these tension 
product operators always can be efficiently converted into a basis operator 
or a sum of two basis operators of the product operator sets Qa and Qb with 
the help of the recursive relation (10), the zero-quantum unitary transforma- 
tion (13), and the single-spin rotation operations: exp(—i9I k ^),fj, = x,y,z. 
As an example, suppose that a product operator Q c is given by 

Qc = 2I lx (g) E 2 0(^3 + h x ) (g) hz ®(\E b - h z ). 
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Then it can be expressed as 

Q c = Ro[(^E 1 + I lz )(g}E 2 (g)(^E 3 + I 3z )(g)E 4 (g}(^E 5 + I 5z ) 

-\e, (g) E 2 (g)(^ 3 + I 3z ) (g) E 4 (g)(^ 5 + h z )]R+ 

where the unitary operator R is determined by the recursive relation (10) 
and the single-spin rotation operations, 

71 7T 

R = exp(i-I ly )exp(i-I 3 y)exp(-iirI 5x ) 

IX 7T 

x exp(-i-I lx ) exp(-inl lz l 4z ) exp(-i-I ly ). 

If the zero-quantum unitary transformation (13) is further used then the 
unitary operator exp(— i6Q c ) can be thoroughly decomposed as a sequence 
of the basic quantum gates: 

exp(-i0Q c ) = R^V 23 (n) + V 35 (n) + C 3 , (e)V 35 (7r)V 23 (7r) 

xV 13 (n) + V 25 (n) + C 2fi (-8/2)V 13 (n)V 75 (ir)R+. 



3. The quantum circuits for the subspace-selective multiple- 
quantum unitary operations 

3.1 The Hermitian diagonal operators 



Two types of Hermitian operators will be used to generate the multiple- 
quantum unitary operators (5). The first type simply consists of the (n + 
1) Hermitian diagonal operators {go, gi, g 2 , ■ g n } for an n— qubit spin 
system. Each of the diagonal operators corresponds one-to-one to one state 
subspace S zq (k). The diagonal operators are generated from the quantum- 
state diagonal operator set {D^}. Their definition is given below. 

go = Do, Do G S zq (0) x S zq (0), 
or in the matrix representation the diagonal operator g is written as 



9o 



Diag(l ,0i, ...,0 



N-l), 
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d(l) 

9l =Z D k ,D k eS zq (l)xS zq (l), 
k 



1 



or g 1 





= Diag(0 o , l h , l Ll , Ll+1 , Ojv_i), 
Z 1 = d(0) = 1, Li = Zi + d(l)-l; 

d(2) 

or 5( 2 = Diag(0 , Oj 2 _i, lj 2 , l La , L2+ i, Ojv-i), 
Z 2 = d(0) + d(l), L 2 = Z 2 + d(2) - 1; 



d(m) 

5- m = D k, D k E S zq {m) x S zg (m), 
k 



or # m 



Diag(0 , ...,0 Z 



-1) J-ir, 



or 



d(0) + + ... + d(m - 1), L m = Z m + d(m) - 1; 
) 

g n = Av-i e 5a, (n) x S^n), 



= Dmg-(0 o ,0i, 







.., 0at_ 2 , 1jv-1 



where again d(m) = (") is dimensional size of the subspace S zq (m); l m 
and L m are diagonal-element indices in the matrix g m ; and S zq (m) x S zq (m) 
is the zero-quantum operator subspace corresponding to the subspace S zq (m) 
whose basis operator set may be simply {|fc)(Z|} with any basis states \k), \l) G 
S zq (m). This zero-quantum operator subspace also includes the diagonal op- 
erator g m . Obviously, the unitary rotation operation G m {6) = exp(—iOg m ) 
(m = 0, 1,2, ...,n) can be expressed as a sequence of the selective rotation 
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operations applied only on the states of the subspace S zq (m): 

G m (9) =H C k (6) (15) 

k=l rrL 

where Ck(9) = C ny k{0) = exp(—i9D k ) given by (11) with the diagonal op- 
erator Dk G S zq (m) x S zq (m) is a selective rotation operation applied only 
to the computational base \tp k ) of the subspace S zq (m), where the computa- 
tional base \(p k ) =<§) {\Ti + afS t ), T\ = |0/> + |1,) and S t = ±(|0,) - |1,» [2]. 
i=i 

It can be seen from (15) that the rotation operation G m (6) of the subspace 
S zq (m) x S zq (m) with m ~ n/2 is an exponential sequence of the selective 
rotation operations with number d(m). But it may be really simplified and 
its efficient quantum circuit may be built up with the elementary propagators 
R klk2 ... km (0) (10) and C^-\9) (12). 

The diagonal matrix g m has d(m) and (2 n — d(m)) diagonal elements 
taking one and zero, respectively. The (2™ — d(m)) zero-diagonal elements 
are divided into two parts by the d(m) one-diagonal elements, and numbers 
for the first part and the second are l rn and (2™ — d(m) — l m ), respectively. 
Note that 2 n = l rn +d(m) + (2 n — d(m) — l m ) is an even number. Then there are 
only two possibilities: (i) all the three numbers l m , d(m), and (2 n — d(m) — l m ) 
are even or (ii) one of the three numbers must be even and the other two 
numbers are odd. For the first case that all the three numbers are even the 
diagonal operator g m can be reduced to the form 

g m = Diag(0, 0; lj m , l Lm ; 0, ...,0) 

= Diag(0, 0; l, m/2 , l Lm/2 ; 0, 0) 2 n-i x2 „-i (g) E n . 

Now the new diagonal operator Diag(0, 0; l im/2 , lh m /2', 0, 0) 2 n-i x2 n-i 
is applied only to the subsystem with the first n—1 qubits of the n— qubit spin 
system. Its dimensional size is 2 n_1 x 2 n_1 , which is denoted in the subscript 
for convenience, instead of dimensional size 2 n x 2™ of the original diagonal 
operator g m . The number of the diagonal elements taking one in the operator 
Diag(0, 0; 1, 1; 0, 0) 2 n-i x2 n-i now is d(m)/2. The unitary operation 
G rn (9) then is simplified (as defined), 

G m (0) = exp[-i0Diag(O, 0; l, m/2 , l Lm/2 ; 0, 0) 2 n-i x2 „-i (g) E n \ 
= exp[-i6Diag(0, 0; l, m/2 , l £m/2 ; 0, 0) 2 n-i x2 n-i]. 



14 



Thus, the unitary operation G m {9) is really applied only to the subsystem 
with the first n — 1 qubits of the n— qubit spin system and hence is reduced 
to a (n — 1)— qubit rotation operation of the (n — 1)— qubit spin subsystem. 
If now the rotation operation G m {9) is still expressed as a sequence of the 
selective rotation operation applied to the (n — 1)— qubit subsystem then 
number of the selective rotation operations, i.e., d(m)/2, in the sequence is 
a half of the original one. Obviously, such a reduction can be further carried 
out. 

It is slightly complicated for another case that only one number is even 
among the three numbers l m , d(m), and (2 n — d{m) — l rn ). In this case there 
are only three possible options: 

(a) l m is an even number. Then d{m) is an odd number. The diagonal 
operator g m can be simplified by 

9m = Diag(0, 0; l Jm/2 , l{ Lm -i)/2] 0, 0) 2 »-i X 2— * ® E n + D Lm 

where the diagonal operator D Lm = Diag(0, 0; lm , Lm -i, U m ; 0, 0), 
and the diagonal operator Diag(0, 0; h m /2, l(L m -i)/2; 0, 0) 2 n-i x2 n-i 
has (d(m) — l)/2 one-diagonal elements and is a (n — 1)— qubit diagonal 
operator. Because all the diagonal operators are commutative to each other 
the corresponding unitary operation G m {6) is decomposed as 

G m {9) = C Lm (6) exp[-i6Diag(0, 0; l Jm/2 , l( Lm _i)/ 2 ; 0, 0) 2 n-i x2 n-i]. 

Therefore, the diagonal unitary operation G m (0) is decomposed into a prod- 
uct of an n— qubit selective rotation operation CL m (6) and a (n — 1)— qubit 
diagonal unitary operation of the n— qubit spin system. 

(b) d{m) is an even number. Then the diagonal operator g m is written as 

9m = Di m + D Lm 

+Diag(0, 0; l(j m+ i)/ 2 , l(L m -i)/25 0, 0) 2 n-i X 2»-i (^) E n . 

The diagonal operator Diag(0, 0; l(/ m +i)/2, -, l(L m -i)/2l 0, 0) 2 n-i x2 n-i is 
a (n — 1)— qubit diagonal operator with (d(m) — 2)/2 one-diagonal elements. 
The corresponding unitary operation G m {6) is expressed as 

G m (0) = C lm (6)C Lm (0) 

x exp[-i0Dia^(O, 0; l(j m +i)/ 2 , l(L m -i)/25 0, 0) 2 n-i X 2»-i]- 
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This shows that the diagonal unitary operation G m (9) now is decomposed 
into a product of two n— qubit selective rotation operations Ci m {9) and 
CL m {0) and a {n — 1)— qubit diagonal unitary operation. 

(c) (2 n — d{m) — l m ) is an even number. The diagonal operator g m is 
simplified by 

9m = D lm + Diag(0, 0; l(« m+ i)/ 2 , U m / 2 ; 0, 0) 2 n-i x2 n-i (g) E n . 

The diagonal operator Diag(0, 0; l(j m +i)/2, — , U m / 2 ; 0, 0) 2 n-i x2 n-i now 
is a (n — 1)— qubit diagonal operator with (d(m) — l)/2 one-diagonal elements. 
The corresponding unitary operation G m {6) can be decomposed as 

G m {9) = C lm (6) exp[-i6Diag(0, 0; l ( z m +i)/2, U m / 2 ; 0, 0) 2 n-i x2 n-i]. 

Thus, the diagonal unitary operation G m (6) now is decomposed into a prod- 
uct of one n— qubit selective rotation operation Ci m (8) and another (n — 
1)— qubit diagonal unitary operation. 

As a summary, in either case the n— qubit diagonal unitary operator 
G m {6) can be reduced to a product of a (n — 1)— qubit diagonal unitary 
operator and two n— qubit selective rotation operators at most. 

The (n — 1)— qubit rotation operation exp[— i9Diag(0, 0; 1 1;0 
, 0) 2 n-i x2 n-i] can be further reduced to the (n — 2)— qubit rotation opera- 
tion exp[— i0Diag(0, 0; 1, 1; 0, 0) 2 n-2 x2 ™-2] which has around d{m)/2 2 
one-diagonal elements, but this reduction may yield extra two (n— 1)— qubit 
selective rotation operations C n -ij m (Q) (the index t m is dependent on l m 
and L m ) at most, so that the unitary operator G m (8) now is a product of a 
(n — 2)— qubit rotation operation and four n— and (n — 1)— qubit selective 
rotation operations at most. This reduction process can continue to the end 
when the diagonal operator g rn is reduced to the final form 

Diag(0, ..,0; 1, 1; 0, .., 0) 2 n^ x2 n- fc (g) E k+1 (g) ... (g) E n (16) 

or 

Diag(0, 0; 1; 0, 0) 2n -, x2 ,- fc (g) E k+l (g) ... (g) E n , (17) 

where k satisfies 2 k w rf(m) and is less than n — 1 because d(m) < 2™ _1 for 
n> 2. The first diagonal operator (16) can form two (n — k)— qubit selective 
rotation operations and the second (17) can generate only one (n — k)— qubit 
selective rotation operation. Since each reduction step can generate two 
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selective rotation operations C^ tm {0) at most the diagonal unitary operator 
G m (8) can be expressed as a sequence of I— qubit selective rotation operations 
Ci t t m {0) (I — k, k + 1, n) with a total number less than 2n. 

The same decomposition procedure as the above can be carried out for a 
general diagonal operator Diag(Oo, Oj_i; 1;, 1l; Ol+i, Ojv-i) that may 
not be in only one zero-quantum operator subspace, e.g., S zq {m) x S zq {m), 
and consequently the diagonal operator may be expressed as a sequence of 
/—qubit selective rotation operations Cj,t m (0) (7 = k, k + 1, n) with num- 
ber less than 2k, where k satisfies 2 k « (L — / + 1) and is less than n. 



3.2 The Hermitian anti-diagonal operators 

Another type of the Hermitian operators used to generate the subspace- 
selective multiple-quantum unitary operators are anti-diagonal Hermitian op- 
erators. They are a generalization of the product operator 2 n ~ 1 I\ x Ii x ---Inx 
[2, 6] and are defined in a matrix form by 



bi 



. 

1 





1 







1 






1 
1 
. 



1 
1 



,6-1 



1 



10 
10 





1 

1 




10 
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o 1 o 



There are 2(2 n -2)+l anti-diagonal operators b k (k = 0, ±1, ±2, ±(2 n -2)) 
in ann— qubit spin system except the two operators |0)(0| and |2 n — l)(2 n — 1|. 
The two operators are anti-diagonal operators and also diagonal operators. 
They are usually assigned to diagonal operators. All these anti-diagonal op- 
erators are symmetrical and Hermitian operators. For every anti-diagonal 
matrix b k all its non-zero elements that take one are located along an anti- 
diagonal line of the matrix. The matrix bo is the main anti-diagonal matrix 
where all its 2 n non-zero elements taking one are located along the main anti- 
diagonal line. The two end points (row, column) of the main anti-diagonal 
line are (0, 2™ — 1) and (2™ — 1, 0) in the matrix bo, respectively. There are 
only (2™ — k) nonzero elements in the matrix b k (or (2 n — 2 > k > 0) 
along the anti-diagonal line of the matrix. Actually, there is a unit difference 
between numbers of matrix element in two nearest anti-diagonal lines in a 
matrix. Thus, nonzero-element number of the anti-diagonal matrix b k+ \ is 
one less than that of the matrix b k , and since the matrix b has 2 n nonzero 
elements the matrix b k has (2 n — k) nonzero elements. The two end points of 
the anti-diagonal line for the matrix b k are (0, 2 n — 1 — k) and (2™ — 1 — k, 0), 
respectively, and for the matrix b- k are (k,2 n — 1) and (2 ra — l,k), respec- 
tively. Denote that x is row coordinate and y column coordinate. Then the 
main anti-diagonal line is given by x = —y + 2 n — 1, and the anti-diagonal 
lines of the matrices b k and b_ k are given by x = —y + 2 n — 1 — k and 
x = — y + 2 n — 1 + k, respectively. Note that an anti-diagonal operator b k is 
symmetric and it has (2 n — k) nonzero (one) elements. If the index k is odd 
then the matrix b k must contain a diagonal element taking one. This diagonal 
element exactly locates at the position (2 n ~ 1 - (k + l)/2,2 n ~ 1 — (k + l)/2) 
along same anti-diagonal line of the matrix b k . However, there is not any 
diagonal element in the anti-diagonal matrix b k with an even index k. 

The unitary operator B k (6) = exp(—i6b k ) built up with any Hermitian 
anti-diagonal operator b k by exponential mapping may be decomposed into a 
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sequence of the basic unitary operations. This may be achieved by express- 
ing the anti-diagonal operator as a sum of basic product operators. Several 
important unitary operators Bk{9) are given below with their explicit decom- 
position. Note that an anti-diagonal operator b±k with k > 2™ _1 first can be 
simplified by 

{b±k)2 n x2 n = (T^ 1 ^ ^ lz ) (gX^i^- 2 ™ -1 )) 2 ™ -1 * 2 ™ -1 

where b±k = {b±k) 2 n x2 n with the subscript 2 n x 2 n indicating dimensional 
size of the matrix bk- Now the index (k — 2 n_1 ) < 2 n_1 and one needs to 
express only the (n — 1)— qubit anti-diagonal operator (&±(fc-2"- 1 ))2"- 1 x2™- 1 in 
the product operator form. Therefore, explicit product operator expressions 
will be given only for those operators b±k with k < 2 n ~ l below. 

(a) The main anti-diagonal operator b . The operator b can be easily 
expressed as 

b = TI lx (g) I 2x (g) ... (g) I nx . (18) 

This is a simple anti-diagonal operator often used in previous papers [2, 3, 5, 
6] . With the help of the single-spin rotation operations and the recursive rela- 
tion (10) the unitary operator B (9) = exp(—i9b ) can be easily decomposed 
into an efficient sequence of one- and two-qubit gates. 

(b) The operator b±\. The product operator expression for the operator 
b\ is slightly complicated. There is a recursive relation for the anti-diagonal 
operator b\\ 

b\ = (&l)2™x2™ = 2™ 1 hxhx---I n -ix 0(2-^™ + hz) 

+ (& 1 ) 2 n-l x2 n-l (g)(\E n - I nz ). 

Using this recursive relation one can express the operator b\ as a sum of n 
commutative product operators: 

h = {\e, + l u ) 0(1^ - I 2z ) (g) ... (g)(i^ n - I nz ) 

+2I lx (g)(^ 2 + hz) (g)(^3 - hz) (g) - ®{\En - Inz) 
+2 2 h x I 2x (g)(^3 + hz) (g)(^4 - hz) (g) - (g)(^n - hz) + 

+2 H - 2 h x (g) hx (g) .. (g) h-2x ®{\En-l + h-lz) (g)(^n - hz) 



19 



+ T~ x h x (g) I 2x (g) ... (g) 7 n _ lx (g)(-£„ + I nz ). (19) 

Since all these product operators in the operator 61 are commutative the 
corresponding unitary operator B\ (9) is decomposed as a sequence of n basic 
unitary operations, 

B 1 (9) = exp[-^(^ 1 + / l2 )(g)(^ 2 -7 2 ,)(g)...(g)(^ n -4,)] 
x exp[-i02/ la (g)(^ 2 + / 22 ) 0(^3 - 4.) (g) - 
(g)(^ n - 4)] 

x exp[-i92 2 I lx I 2x (g)(^ 3 + 4.) (g)(^4 - 4.) (g) - 
®{\E n - I nz )\ x 

x exp[-t9T- 2 I lx I 2x ...I n . 2x (g)(^„-i + /„_u) (g)(^n - 4,)] 
x expH^"" 1 /^^. ..I n _ lx (g)(^„ + J n2 )]. (20) 

In a similar way, the anti-diagonal operator 6_x also can be expressed as a 
sum of n commutative product operators, 

6_i = (l^i - 7 l2 ) <g)(i£ 2 + J 22 ) ... ®{\E n + J n2 ) 

+2J lie <g>(±£ 2 - I 2z ) (g)(§£ 3 + / 32 ) <S> - ®{\En + 4.) 

+2 2 hJ 2x (8)(^ 3 - 4.) ®(^4 + 4z) ® - (g)(|^n + hz) + 

+2"~ 2 / 1:c (g) 7 2x (g) ... (g) J n „ 2:c (g)(i£? n _i - 4-1.) (g)(|K + 4.) 
+2"" 1 /!, (g) 7 2a; (g) ... (g) / n _ lie <g)(±£ n - I nz ), 

and corresponding unitary operator B_i(9) therefore is decomposed as a 
sequence of n basic unitary operations. 

(c) The operator b±k with k — 2 l (I — 1,2,..., n— 1). First consider 
the operator b k with an even index k. The number of nonzero (one) ele- 
ments of the matrix b k along the anti-diagonal line is 2 n — k. If the in- 
dex k is even then so is 2 n — k. If now the matrix bk is blocked by a 
2x2 submatrix one can see this blocked matrix bk is still an anti-diagonal 

blocked matrix, and the nonzero blocked submatrix is 2I X = 

along the anti-diagonal line. Therefore, the matrix 6 fc may be written as 
(bk) 2 n x2 n — (&fc')2 n - 1 x2™- 1 (g)27 n:r with index k' = k/2, and by is also an anti- 
diagonal matrix. If k' is still even then the matrix bk can be further written 
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as (bk)2 n x2 n = (&fc")2"- 2 x2"- 2 (£) 1I n -\ x 2I nx with index fc" = fc/4. Gener- 
ally, for the operator bk with k — 2 l (I — 1,2, — 1) number of nonzero 
elements is (2 n — 2 l ) on the ant i- diagonal line. Then the operator b 2 i can be 
reduced to the form 

b 2 i = (6i) 2"- ; x2 n -' (/0 2 l In-l+lxIn-l+2x---Inx- 

(21) 



Particularly the matrix 62 can be written as 62 = (^i)2 n - 1 x2 n - 1 (£)2I n x- Be- 
cause the operator (6i) 2 n-!x2™- i can be further expressed as a sum of (n — I) 
commutative product operators, as shown in (b), the operator b 2 i now is 
written as a sum of (n — I) commutable product operators and hence the 
unitary operator B^{9) = exp(— iOb^) with k — 2 l (I — 1, 2, n — 1) can be 
decomposed as a sequence of (n — I) basic unitary operations. In a similar 
way, the anti-diagonal operator b-k with k = 2 l can be reduced to the form 



b_ 2 i — (b-i) 2 n-i x2 n-i (x) 2 l I n _i + i x I n _i +2x ...L, 



1 nx j 



where {b^\) 2 n,-i x2n -i can be further expressed as a sum of (n — I) commutative 
product operators. Therefore, the unitary operator = exp(— i9b-k) 

with k = 2 l can also be decomposed as a sequence of (n — I) basic unitary 
operations. 

Generally, an anti-diagonal operator bk with an even index k = 2 kl ~ 1 + 
2 hl - 2 + ... + 2 hl can be simplified by 

bk — (bk') 2 n-k lx2 n-k 1 (X) 2 kl I n _k 1+ i x I n _k 1+2x ... I nx (22) 



with the odd index k' = 2 k '-^ kl + 2 kl ~ 2 - kl + ... + 2 k ' 2 - kl + 1. 

(d) A general operator 6 fc . In (c) it is shown that an anti-diagonal oper- 
ator b k with an even index k can be reduced to another lower- dimensional 
anti-diagonal operator with an odd index. Here consider the operator b k with 
an odd index k. The operator bk always can be written in the form 

(frfci/2)2"- 2 x2™- 2 <8> 2ira-lx ®{\E n + Inz) 

n . ) +{bi 1 ) 2 n-i x2 n-i ®(\E n - I nz ), if ki is even, , . 

1 " ,2 " x2 " " (6i 1 /2)2- a x2- a ®2/ B _ be ®(iS B -/„) U,) 
+(6fc 1 ) 2 ™-ix2»-i (8)(i^n + 4*), if h is even, 

where ki = (k — l)/2 and Zi = (A; + l)/2. The number of nonzero elements 
keep unchanged before and after the reduction: (2 n — k) = (2 <n ~ 1 ) — hi) + 
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(2( n x ) — The n— qubit anti-diagonal operator (bk)2 n x2 n now consists of 
one term containing (n — 1)— qubit anti-diagonal operator (6^) 
(^fci)2™- 1 x2 n - 1 ) and another term containing (n — 2)— qubit anti-diagonal op- 
erator (&fe 1 /2)2"- 2 x2"- 2 ( or (^i/2)2"- 2 x2"- 2 )- Note that the two terms are com- 
mutable to each other. This is an advantage of the decomposition based on 
the recursive relation (23). Image that using (23) an n— qubit anti-diagonal 
operator at the first step is decomposed as a sum of two (n — 1)— qubit 
anti-diagonal operator terms, then at the second step the two (n — 1)— qubit 
anti-diagonal operators are decomposed as four (n — 2)— qubit anti-diagonal 
operator terms, and so on, in the final the n— qubit operator bk would be a 
sum of 2 n_1 commutative basic product operators if the decomposition could 
be carried out to the final (n — l)th step. However, the recursive relation (23) 
shows that although the n— qubit operator bk is first decomposed as a sum of 
two (n — 1)— qubit anti-diagonal operator terms, one of the two (n — 1)— qubit 
anti-diagonal operators can be further reduced to only one (n — 2)— qubit 
anti-diagonal operator term instead of two terms. Then term number of the 
basic product operators in the completely decomposed operator 6 fc is not re- 
ally more than 2 n ~ l . For some specific cases one may use conveniently the 
recursive relation (23) to decompose an anti-diagonal operator bk as a sum 
of polynomially many basic product operators. Take the anti-diagonal op- 
erators b± k with index k = T ± 2 m (n > r > m — 0, 1, 2, ...,n — 1) as an 
example. 

For the index k = 2 r ± 2 m the operator bk is first reduced to the form 

(bk)2 n x2 n — (&2 r - ™±l)2™- m x2 n - m ®2 m In-m+lxIn-m+2x---Inx, &S shown in (22). 

Then consider the anti-diagonal operator (&2 ! +i)2*x2* with t = n — m and 
I — r — m. It follows from (23) that there is a recursive relation for the 
operator (b 2 i +1 ) 2 tx2 t - 

(^2 i +l)2*x2 i = (^2'- 1 +l)2*- 1 x2*- 1 (^Xt^ ~~ 

+(b 1 ) 2 t-i X 2t-i 2 l ~ 1 I t _ l+lx I t _i +2x ---It-ix (^Xt^ + 

This is because in (23) k\ = (k — l)/2 = 2 l ~ 1 and {b 2 i-i / 2 ) 2 t-2 x2 t- 2 can be 
further reduced to the operator (bi) 2 t-i shown in (c). This relation 

directly leads to the product operator expression for the operator (b 2 i +1 ) 2 t x2 t : 

(&2 i +l)2'x2 t = (&l)2'- i x2'- i ® ^ l ~ 1 h-l+lxh-l+2x--h-lx <S>( 2 E t + hz) 
+ (b 1 ) 2 t-l x2 t-l <S>2~ It-l+lxh-l+2x---h-2x (g>( 

®C 2 E t - I tz ) 
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+ (h) 2 t-l x2 t-l<S}2 1 Z h-l+lxh ®(§£ t _ 2 + / t _ 2z ) 

(8)(|^_i - i t -u) <g>(lE t - I tz ) + 

-I ® 2It-l+lx ®(\Et-l+2 + h-l+2 Z ) ®{\E t 

-1+3 — h-l+3z) 

<g)...<g>(lE t -I tz ) 

+(h) 2 t-i x2 t-i (g)(§£ t _ m + h-i+iz) ®{\E t - l+ 2 - I t -i+2z) <S> ••• 
®{\E t -I tz ) 

+ (^2)2'- i x2'- i &){\E t -l+l - h-l+lz) (£){\E t -l + 2 - h-l+2z) <S> ■■■ 

<S)(lE t -I tz ). 

The operator (& 2 i +i)2 t x2 t now is a sum of (/ + 1) commutative operators. 
Since the operator ipi) 2 t-i x2 t-i can be expressed as a sum of {t — I) com- 
mutative product operators, as shown in (b), and (b 2 ) 2 t-i x2 t-i as a sum °f 
(t — l — 1) commutative product operators, as shown in (c), then the operator 
(b 2 i + i) 2 t x2 t is a sum of (I + l)(t — I) — 1 commutative product operators. 
Now taking t = n — m and I — r — m one sees the operator (&fc) 2 ™ x2 ™ with 
k = 2 r + 2 m is a sum of (r — m + l)(n — r) — 1 commutative product oper- 
ators. Thus, the unitary operator B^iO) can be decomposed as a sequence 
of (r — m + l)(n — r) — 1 basic unitary operations. In an analogous way, it 
can be shown that the unitary operator can be also decomposed as 

a sequence of (r — m + \){n — r) — 1 basic unitary operations. 

From (23) there is also a recursive relation for the operator b 2 i_ x : 

(&2'-l)2*x2* = (^-i-l)2'- 1 x2*- 1 ^)(^E t + I tz ) + (b 2 l-i) 2 t-i x2 t-i (^{^Et- I tz ) . 

Then using this relation the product operator expression for the operator 
(&2*-i)2*x2* is given by 




(g)(^-, + 3 + It-l+3z) - + /**)■ 
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Because the operator (b 2 i-m) 2 t-m x2 t-m ( m = 1,2,...,/ — 1) can be expressed 
as a sum of (t — I) commutative product operators the operator (&2 ; -i)2*x2* 
is clearly a sum of l(t — I) + 1 commutative product operators. Then the 
operator b k with k = T — 2 m is a sum of (r — m){n — r) + 1 commutative 
product operators and hence its unitary operator B k {6) can be decomposed 
as a sequence of (r — m){n — r) + 1 basic unitary operations. 

The above results suffice to construct an efficient subspace- selective 
multiple-quantum unitary operator (5) that transfers completely an unknown 
state from any subspace S zq {m) {m ^ n/2) to the largest subspace S zq {n/2) 
of the Hilbert space of an n— qubit spin system. 

A general anti-diagonal operator also can be decomposed using the recur- 
sive relation (23) in an analogue way. But by using only the recursive relation 
(23) it is usually not convenient to obtain an efficient decomposition for the 
unitary operator B k {6) = exp(—i9b k ) with a general anti-diagonal operator 
bk- It can be proved that a general anti-diagonal operator can be converted 
into another anti-diagonal operator with a different index by a proper uni- 
tary transformation, and for a general anti-diagonal operator bk with an odd 
index k = 2 kl - 1 + 2 kl - 2 + ... + 2 kl + 1 (n - 1 > fc,_! > k_ 2 > ••• > h > 1) its 
unitary operator Bk{9) = exp(—i9bk) can be generally expressed as 

B k (e) = U k exp(-iOb kl )U+ (24) 

where the operator b k \ is a symmetric and Hermitian anti-diagonal operator 
similar to the anti-diagonal operator bi and U k is a unitary operator depen- 
dent on the index k. The anti-diagonal line of the operator b kl is the same 
as that of the operator bi, but number of nonzero elements taking one along 
the anti-diagonal line in the operator b kl is only (2™ — k) instead of (2™ — 1) 
of the operator b\. The (2 n — k) nonzero elements locate symmetrically at the 
center of the anti-diagonal line and each of two ends of the anti-diagonal line 
has (k — l)/2 zero elements. Note that the index k is odd. Just like b± the 
symmetric matrix b ki has a diagonal element at position (2 n_1 — 1, 2 n_1 — 1) 
and hence the operator b kl also contains the diagonal operator D 2 ^-i^i. The 
diagonal operator D 2"- 1 — i is commutable with both the operators b\ and b k \. 
The unitary operator B k i(0) = exp(—i9b k i) can be decomposed efficiently. 
Just like the subspace-selective multiple-quantum operator (27) (see next 
section) the Hermitian operator b k i can be expressed as 

(Pkl - ^2»-i-l) = ( b l ~ £>2»-l-l)] + 

= 9k( b i - D 2 n-i_ 1 ) + (6i - D 2 n-i_ 1 )g k 
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where the diagonal operator g k is given by 

g k = Diag(0 , 0( fc _i)/ 2 _i; l( fe _i)/ 2 , l( 2 ™-i_i)_i; 0( 2 ™-i_i), Ojv-i)- 

As shown in section 3.1, the unitary diagonal operator G k {0) = exp(—i9g k ) 
can be decomposed efficiently into a sequence of the I— qubit selective rotation 
operations C^ tm (0) (I = 1,2,..., n) with number less than 2n. Since both 
the two operators (b~ki — D 2 n-i_i) and (b\ — D 2 ™-i-i) do not contain any 
diagonal operator components and both the two operators g k {b\ — D 2 n-i_ 1 ) 
and (&i — D 2 n-i_ 1 )g k have not any overlapping nonzero matrix element it 
follows from the unitary transformation (38) and the decomposition formula 
(41) in next section that the unitary operator Bk\{0) = exp(—i9bki) can be 
efficiently decomposed as 

B kl (6) = C 2 n-i_ 1 (e)exp[-iHbki-D 2 n-i_ 1 )} 

= C 2 „-i_ 1 (0)(ex P H^(6i - D^-^/LjGkin) 

x exp[i^(6i - L» 2 n-i_ 1 )/L]G fe (7r)- 1 ) L + O^ 1 ) 

= C 2 „-i-i(0)[ex P (-i^6 1 /L)G fc (7r) 

x exp(i^6 1 /L)G jfc (7r)- 1 ] i + O^- 1 ). (25) 

For a modest number L = e -1 this expansion of the unitary operator B kl (9) 
can fast converge. As shown in (b), the unitary operator exp(— i\Qb\jV) is a 
sequence of n basic unitary operations. Then the unitary operator B k i{6) can 
be decomposed as a sequence of the basic unitary operations with number 
less than 6e~ 1 n. 

It can be proved that the unitary operator U k in (24) with the index 
k = 2 kl -i + 2 fc '- 2 + ... + 2 fcl + 1 (n - 1 > ki-! > k^ 2 > ... > h > 1 and 
1 < I < n — 1) can be written as 

{exp(if b h ) exp(if &_ ja ) exp(if b h ) 
x... x exp(if b-j^), if I - 1 is even, 
exp(i!& jl )exp(if&_j a )exp(if& J 3) 
x... x exp(i|6 J; i ) exp(i|& ), if / — 1 is odd, 

where the indices j\ = 2 kl - 1 ~ 1 , j 2 = 2 kl - 2 ~ 1 , = 2 fel_1 . Since the unitary 

operator exp(— i6b±j) with j — 2 l (I — 1, 2, n — 1) can be decomposed as a 
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sequence of (n — l) basic unitary operations, as shown in (c), then the unitary 
operator exp(i^b±j m ) with index j m = 2 k '- m ~ 1 in (26) can be decomposed 
into a sequence of (n — ki- m + 1) basic unitary operations. Therefore, number 
of the basic unitary operations in the unitary operator Uk (26) is (n — ki-i + 
1) + (n — ki-2 + 1) + ••• + (n — ki + 1) if I — 1 is even or (n — ki-i + 1) + (n — 
h-2 + 1) + ■•• + (n - h + 1) + 1 if I - 1 is odd. Note that (n - fcj_i + 1) + (n - 
ki_2 + 1) + ... + (n — /ci + 1) + 1 < (I — l)n + 1 < n 2 . The unitary operator Uk 
(26) then can be decomposed into a sequence of the basic unitary operations 
with number less than n 2 . 

Therefore, the expansion (25) of the unitary operator exp(— iObki) and 
the decomposition (26) of the unitary operator Uk show that the unitary 
operator B k (9) (24) built up with a general anti-diagonal operator b k can 
be expressed as a sequence of the basic unitary operations with number less 
than 2n 2 + 6e~ 1 n, and quantum-circuit complexity for the unitary operator 
is 0(2n 2 + Qe~ 1 n). 

3.3 The subspace-selective multiple-quantum operators 

The Hermitian multiple-quantum operator Q pm in the subspace-selective 
multiple-quantum unitary operator U pm (9) = exp(—i9Q pm ) (5) can be gen- 
erated by the anti-commutator of the diagonal and anti-diagonal operators: 

d(m)-l 

Qpm = ^ Qpml = - [bk,9m]+ = -(Ombk + b k 9m), (27) 
1=0 

where the anti-diagonal operator bk needs to be chosen properly and the 
diagonal operator g m e S zq (m) x S zq (m) so that the p— quantum unitary 
operator U pm (9) built up with the Hermitian p— quantum operator Q pm is 
selectively applied on both the state subspace S zq {m) and another larger 
subspace S zq {m + p). It needs first to show how to choose the anti-diagonal 
operator bk to generate the Hamiterian operator Q pm . Here consider only the 
case < m < n/2. For the case n > m > n/2 the multiple-quantum operator 
Qpm can be constructed with the operator in place of the operator bk 
and the final result is similar. For convenience the anti-diagonal operator 
b k = {bk)NxN (N = 2 n ) is written as 

bk = Adiag(l[o t N-k-i], l[i,JV-fc-2], 1[jv— fc— 1,0] j 0[jv-fc,jv-i], ■■■■> 0[jv-i,jv-fe])- 
Suppose that the operator 6 fc is chosen properly so that the operator g m bk is 
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given by 

g m b k = Diag(0 , O im _i; l Jm , l Lm ; Lm+ i, Ojv_i) 
xAdiag(l[ 0j N-k-i], l[i,JV-fc-2], •••> l[N-k-i,o}] 

0[jV-fc,AT-l], 0[AT-l,JV-fe]) 

= Adia^(0[ 0j jv-fe-i], •••,0[ Jm _i j jv-fe-i-z m +i]; 

l[Z m ,JV-fc-l-J m ], •••) l[L m ,JV-fc-l-L m ]; 

0[L m +l,AT-fe-l-t m -l], •••> O[jv-fc-l,0]! 0[JV-fc,JV-l], 0[7v__i 5 AT_ fc ]), 

and the operator is written as 

bkdm = Adiag(l[ 0j N-k-i], l[i,w-fe-2], •••> l[Jv-fc-i,o]; 0[jv-fc,iv-i], 0[jv-i,jv-fc]) 
xDiag{0 , Jm _i; l, m , l Lm ; Lm +i, Otv-i) 
= Aiia#(0[ ,7v-fc-i], 0[Ar-fc-i-L m -i,L m -i]; 

l[JV-fc-l-L m ,L m ], l[JV-fc-l-J m ,im]; 

0[JV-fc-l-U+l,im-l]> 0[JV-fc-l,0]; 0[JV-fc,JV-l], •••) 0[jv-l,jV-fc])- 

Obviously, the index TV — 1 — k must be greater than the index L m . Both 
the operators g m bk and b k g m are also anti-diagonal operators. Their nonzero 
elements taking one are also along the same anti-diagonal line of the matrix 
bt- Each of the two operators has only d(m) nonzero matrix elements, which 
number is exactly dimensional size of the state subspace S zq (m). The opera- 
tor (g m bk + bkgm) is clearly a symmetric and Hermitian operator and can be 
written as 

(g m bk + b k g m ) 
= Adia^(0[ 0j iv-fe-i], 0[z m _i )JV -fc-i-z m+ i]; 

l[l m ,N-k-l-l m ]iQ[l m +l,N-k-l-l m -i\, ■■■,Q[N-k-l-l m -l,l m +i\'> 
l[N-k-l-l m ,l m ]'iQ[N-k-l-l m +l,l m -Hi ■■■■>Q[N-k-l,0]', 
0[JV-fc,JV-l], •••) 0[jv-l,JV-fc]) 

+Adiag(0[ 0jN - k - 1 ], 0[i m ,N-k-i-i m ], 

l[J m +l,JV-fc-l-im-l]) Q[l m +2,N-k-l-l m -2], ■■; Q[N-k-l-lm,-2,l m +2], 

l[N-k-l-l m -l,l m +l],Q[N-k-l-l m ,l m ]i ■■■■> 0[JV-fc-l,0]; 
0[JV-fc,JV-l], •••) 0[AT-l,iV-fc]) + 

+AZia0(O[ O ,jv-fc-i], 0[L m -i,Ar_fc-i-L m +i], 
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,N-k-l-L m ]j^[L m +l,N-k-l-L m -l]7 ■■■■> Q[N-k-l-L m -l,L.. 
k-l-L m ,L m ]',Q[N-k-l-L m +l,L m -i\, 0[JV-fc-l,O] ! 
k,N-l], 0[N-l,N-k])- 



+ 1]' 



Therefore, the operator (g m bk + bug™) can be expressed as a sum of d{m) 
commutative and symmetric anti-diagonal operators at most, each of which 
has only two elements taking one along the same anti-diagonal line of the 
matrix bj,. Such an anti-diagonal operator can be written in terms of the 
usual computational basis {|<^ fe )} : 

Adiag(0[ O! N-k-i], Q[i m +i'-i,N-k-i-i m -i'+i\i 

l[l m +l',N-k-l-l m -l'],Q[l m +l'+l,N-k-l-l m -l'-i\, Q[N-k-l-l m -l'-l,l m +l'+i\i 

l[N-k-l-l m -l',l m +l']l Q\N-k-l-l m -V+\,l m +l'-l], — 1 0[JV-fc-l,0] ! 
0[N-k,N-l], ■■; 0[JV_l,JV-fc]) 
= IPlm+l') {PN-k-l-lm-l'l + l^N-h-l-lm-l') ( L Pl m +l'\l 

where the index I' = 0, 1, ...,L m — l rn (L m — l m — d(m) — 1) and the base 
\<Pi m +i') belongs to the state subspace S zq {m). If both the operators g m bk an d 
bkg m have common nonzero matrix elements then the operator (g m bk + bkg m ) 
is a sum of commutative and Hermitian anti-diagonal operators with number 
less than d(m). If the basis state \ip N _ k _ 1 _i m _ii) belongs to another subspace 
S(M m +p) this anti-diagonal operator is a state-selective p— quantum opera- 
tor that applies only on both the two basis states \ip lm+l ,) and \(PN-k-i-i m -i')- 
Now it is required that both the operators g m bk an d b k g m have not any com- 
mon nonzero matrix elements so that the operator {g m bk + bkg m ) is a sum of 
the commutative anti-diagonal operators with number d(m) exactly. Then 
the index N — 1 — k > 2L m , and the operator (g m bk + bkgm) can be written 
in terms of the usual computational basis 



Note that all the d(m) basis states \<Pi m+ i>) (I' = 0, l,...,d(m) — 1) in (28) 
belong to the subspace S zq (m). If now all the d(m) basis states |y%-fc-i-; m -2') 
in (28) belong to the subspace S zq {m + p) whose dimensional size is larger 



-(g m bk + b k g m ) 

d(m)-l 



'N-k -l-l m -l'\ ' WN-k-l-lm-l ■)(Vi rn+ i>\)- (28) 



l'=0 
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than that of S zq (m), the operator ^(g m bk + bkg m ) is a subspace-selective 
multiple-quantum operator which applies only to both the two subspaces 
Szq{ m ) and S zq (m + p). This can be seen more clearly by comparing the 
operator (28) with the operator (27) and the multiple-quantum operators 
(3) and (5). Therefore, the condition that the operator \{g m bk + b k g m ) is the 
subspace-selective p— quantum operator Q pm (5) selectively applied to both 
the subspaces S zq (m) and S zq (m + p) is that in addition to iV — 1 — k > 2L m 
the index k for the anti-diagonal operator bk must satisfy 

l m+P <N-k-l-l m -l' < L m+P , I' = 0, 1, .., d(m) - 1. (29) 

Note that L m = l m + d{m) — 1 and L m < l m+p for < m+p < n/2 and p > 1. 
The condition N — 1 — k > 2L m always holds if the condition L m < l m+p < 
N — k — 1 — l m — (d(m) — 1) holds. Thus, the condition (29) is a general 
condition to determine the proper index k for the anti-diagonal operator bk 
that is used to construct the subspace-selective p— quantum operator Q pm 
(27). 

Consider the n— qubit spin system with an even qubit number n. Then 
the largest subspace for the system is S zq (n/2). Let S zg (m + p) = S zq (n/2). 
Then the condition (29) is rewritten as 

/n/2 <N-k-l-l m -l' < L n/2 , I' = 0, 1, .., d(m) - 1. (30) 

The condition (30) shows that l n /2 < N — k — 1 — l m < L n / 2 if V = and if 
/' = d{m) — 1 then Z n / 2 < N — k — 1 — L m < L n / 2 , and since N — k — \ — l m > 
N- k-1- (l m + 1) > ... > N- k- 1 -L m one has for I' = 0, 1, ...,d(m) - 1, 

L/2 <N-k-l-L m <...<N-l-k-l m < L n/2 . 

Therefore, the lower bound (/c m ) m in and upper bound (k m ) max for the index 
k = k m of the subspace S zq (m) for m — 0, 1, n/2 — 1 are given by 

(^m)min = N — l m — L n / 2 — 1 and (/c m ) max = N — L m — l n / 2 — 1, 

and hence the index k m is bounded on by 

N -l m - l n / 2 - d(n/2) < k m < N - l m - l n/2 - d(m). (31) 

The condition (31) shows that range of the index k m is equal to Ak m = 
(^m)max _ (^m)min = d(n/2) — d(m) . For a different state subspace S zq (m) 
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the distance Ak m is different, and the maximum and minimum Ak m are 
AA; = d(n/2) - d(0) = ( ™/ 2 ) ~ 1 &nd AA: ™/ 2 - 1 = d< ^ n / 2 ^ ~~ d ^ n / 2 ~ l ) = 
™+2 iy ) ' res P ec ^ ve ^' anc ^ moreover AA; > Aki > ... > Ak n / 2 -i. Note 

that 2l n/2 + d(n/2) = £™ =0 d ( m ) = E m=0 ( ^ ) = 2 ™ = ^ The condition 
(31) is reduced to the form 

L/2 -l m <k m < l n/2 - l m + d(n/2) - d(m) (32) 

where l m = d(0) + d(l) + ... + d(m — 1) and Iq = 0, and l n / 2 > l n /2-i > 
... > li > Iq. Now using the condition (32) one can determine the index 
k m for the desired operator bk m . Suppose that the index k m = 2 n ~ 1 . Since 
N/2 = l n / 2 + d(n/2)/2 > l n / 2 the first inequality in the condition (32) always 
holds: k m = N/2 > l n / 2 — l m . The second inequality is reduced to the form 

d(n/2)/2 -l m - d{m) > 0. (33) 

Therefore, the operator bk m with k m = 2 n_1 can be used to construct those 
p— quantum operators Q pm with index m — 0,1, ...,m where the maximum 
index m is determined from the inequality (33): 

1=0 

As shown in section 3.2, the anti-diagonal operator b km with index k m = 2 n ~ 1 
can be explicitly expressed as 

h m = [\e, + I lg ) (g)2 n - 1 I 2x I 3x ...I nx 

and the corresponding unitary operator B km {9) therefore is a basic unitary 
operation, 

B km {6) = exp[-t0(^E 1 + I lz )(g)2 n - 1 I 2x I 3x ...I nx ]. (34) 

This quantum unitary operator is used to build up a subspace-selective 
p— quantum unitary operator U pm {6) (5) that can transfer any unknown state 
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that is in one of the subspaces S zq (0), S zq (l),..., S zq (m Q ) into the largest sub- 
space S zq {n/2). 

Next consider the situation m < m < n/2 — 1. Since l n / 2 — l m — d(m) + 
dim + 1) + ... + d(n/2 — 1) the condition (32) is rewritten as 

d(m)+d(m+l)+...+d(n/2-l) < k m < d(m+l)+d(m+2)+...+d(n/2). (35) 

In particular, for m = n/2 — 1 the condition (35) is written as 



n/2 - 1 ) " K ' 2 - x ~ ( n/2 



< fcn/2-1 < ( m/0 ). (36) 



By the minimum distance Ak n / 2 -i = ( ^2 1 one ma ^ °^^ am ^ ne m dex 

fc m that satisfies (35) for each subspace S zq {m) with m < m < n/2 — 1. The 
minimum distance A/c n / 2 -i is approximated by using the Starling's formula 
n! w y/2im{n/e) n for a large n, 



AA;n/2 ~ 1 - nT2 [n/2 ) * nT^"^ 
and 

log 2 AA; n/2 _i « n - toga^n-v/ra^^l + ^)). 



Denote n = \}og 2 {\n^Jn^J^{l + |))J as integer part of log 2 (|n v /ny / |(l + 
£)). Then 2 n - n °- 1 < Ak n/2 -i < 2 n - n ° and 2™° ~ n 3/2 . The minimum index 
(^n/2-i)min that satisfies (36) is approximated by the Starling formula 



(fcn/2-l)min - ( n /2-} ) ~ \l ~~~T~n~nz^ n 



n 1 



7r n + z vn 



and 



1 °g 2 (" /2 _ 1 )«^- 1 0g2(v^y|(l + ^)- 



Let A; = ^(v^a/I ■ Then 2 ™" fc °" 1 ^ (^/2-i)min < 2"" fe " and 
2 fe( » ~ n 1//2 . Although taking fc n / 2 -i = 2 n ~ k ° the first inequality in (36) can be 
satisfied, the second inequality in (36) could not be satisfied. Now suppose 
that the index fc n /2-i = 2 n ~ k °~ 1 + /j,2 n ~ n °~ 1 , u is an integer to be deter- 
mined. Then there is an integer u such that the index fc n / 2 -i satisfies the 
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condition (36). There is always an integer fx such that 2 n k ° 1 + fi 2 n n ° 1 < 
(fc»/2-i)min ^ 2"- fc «- 1 + (/i + l)2 n ""o- 1 . The integer /i equals L((fcn/2-i)min- 
2 n - fc °- 1 )/2 n - n °- 1 J. Because 2 n ~ k °~ 1 + /i 2 n - no " 1 < (fc n / 2 -i)min and 2 n ~ n °- 1 < 
A/c n / 2 _i there must be 2"- fc °" 1 + (/x + l)2 n ~ n ^ 1 < (k n/2 ^) min + Ak n/2 ^ = 
d(n/2). By taking /i = /i + 1 the index fc n/2 _i = 2 n - fco ~ 1 + fi2 n ~ no ' 1 sat- 
isfies the condition (36). Generally, suppose that there is an integer /i 
such that 2 n ~ k °- 1 + ^2 n ~ n °- 1 < d(m) + d(m + 1) + ... + d(n/2 - 1) < 
2n _ fco _i + ^ + 1 ^ ra _ no _i_ Because 2 n-n -i < AA; n/2 -i < A/c m there holds 

2 r i -feo-i + ( /io + 1 ) 2 ™"™o- 1 < d(m) + d(m+l) + ... + d(n/2-l)+Ak m = d(m + 
1) + d(m + 1) + ... + d(n/2). Therefore, the index k m = 2 n - k °~ 1 + / u m 2 n " n °- 1 
with /j, m — pL + 1 satisfies the condition (35). For m < m < n/2 — 1 the 
upper bound for the index k m is obtained from the condition (35) by taking 
m = m + 1 : 

k m < d(m + 2) + d(m + 3) + ... + d{n/2) 

= l n/2 + d{n/2) - d(0) - d(l) - ... - d(m ) - d(m + 1) 
< N/2. 

mo+l mo 

This is because ^ d(l) > \d(n/2) and Yl d(l) < |d(n/2), as shown in 

i=o ' 1=0 

(33). Since k m < 2 n ~ 1 the index k m can be expanded as a binary number: 

h. _ on-fco-l _i_ on-no-l 

= a n . 2 2 n - 2 + a n _ 3 2 n - 3 + ... + a n . no . 1 2 n - n °- 1 (37) 

where a\ = 0, 1 (/ = n — n — 1, n — n , n — 2). The anti-diagonal operator 
6 fem with the index (37) then can be expressed as, as can be seen in (22), 

(bk m )2 n x2 n — (bk' m ) 2 " +i x2 n +i (^) 2™ n ° 1 I no +2xIn +3x---Inx 

where k' m = a„_ 2 2 n °~ 1 + a„_ 3 2 n °~ 2 + ... + a n _ no _i2° and the anti-diagonal 
operator by has a dimensional size 2 no+1 x 2 n ° +1 . Because 2™° ~ n 3//2 the 
operator (&*.' ) 2riQ + l x 2 n 0"'" 1 Colli be expressed as a sum of ~ n ' commutative 
product operators at most using the recursive relation (23) in section 3.2. 
Therefore, the operator bk m is also a sum of ~ n 3 ^ 2 commutative product 
operators at most. 

Furthermore, if using the general decomposition (24) for a general anti- 
diagonal operator the unitary operator exp(— i0b km ) with the index k m (37) 
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can be decomposed into a sequence of the basic unitary operations with a 
complexity only O(2(n + l) 2 + 6e _1 (n + 1)), which is ~ O (login). 

For an n— qubit spin system with an odd qubit number the largest sub- 
spaces are S zq {n/2 — 1/2) and S zq {n/2 + 1/2). Any state in one of the 
two largest subspaces can be converted unitarily into another largest sub- 
space by a unitary transformation with n single-spin unitary operations: 
U = YYk=i ex P( — iwlkx)- This unitary transformation is also available for the 
complete quantum-state transfer between any pair of symmetric subspaces 
S zq {m) and S zq (n — m) in a general n— qubit spin system. Thus, it needs 
only to consider those subspaces S zq (k m ) with < k m < (n — l)/2. Take 
S zq {m + p) — S zq (n/2 — 1/2). Now the minimum distance A/c( n _i)/ 2 -i = 

d(n/2 — 1/2) — d((n — l)/2 — 1) = ^ ^ l)/2 ) an< ^ * S a PP rox ^ ma ^ ec ^ 
by 



log 2 Afc(„_i)/2-i « n-log 2 {J^— 



3 (1 +n) v / ^T 



4 n 

x v (i+ ^W [i+ ^r ](n_i)} - 

Now let n = Llog 2 vff ^ ^ + ^V^t 1 + Sl^J- Note that 
(l + -) n ~ e for a large integer n. Then 2™° ~ n 3 ^ 2 . Therefore, for the case that 
the qubit number n is odd the operator bk m with k rn = 2 n ~ 1 is still used to 
construct those p— quantum operators Q pm (27) with index m — 0, 1, ...,m , 

mo 

but now the maximum index m is determined from the inequality: Yl d(l) < 

1=0 

\d((n — l)/2), and when m < m < n/2 — 1/2 the index k m is still given by 
(37). The complexity of quantum circuit of the unitary operator exp(— i0b km ) 
is also ~ 0(log 2 n). 

For a general subspace S zq (m + p) instead of the largest subspace in 
an n— qubit spin system the index k of the operator b k used to construct 
the p— quantum operator Q pm (27) is generally determined from the general 
condition (29), and the unitary operator exp(—i9bk) is decomposed as a se- 
quence of the basic unitary operations by the general decomposition (24), 
and the complexity of quantum circuit of the unitary operator exp(— i9bk) is 
0(2n 2 + 6ne- r ). 

Once the unitary operator exp(— i6b km ) with the anti-diagonal operator 
6 fcm is expressed as an efficient sequence of the basic unitary operations the 
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subspace-selective p— quantum unitary operator U pm (9) (5) can be easily de- 
composed as an efficient sequence of the basic unitary operations. A gen- 
eral unitary transformation with a sequence of selective rotation operations 
{Cfc(^fc)} can help simplify the decomposition. The unitary transformation 
has been given in Ref . [6] : 

U o (0 , #1, 9 m -i)p I (t)U o (9 , 0i, m _i) 1 

m— 1 m— 1 

= Pi(t) ~ [pi(t),J2 (l-co&9 k )D k ] + + i\p I (t),^D k sme k ] 

k=0 k=0 

m-lm-1 

+ [( X ~ cos6 k)( l - c °s0/) + sin0 fc sin0 / ] J D fe p / (t)A 

k=0 1=0 

TO— 1 

+ ^[sm0 fc (l-cos0 z )-sin0 / (l-cos0 fc )]( J DfeP/(t)A-AP/(t) J Dfc) (38) 

l>k=0 

m' — l 

where the diagonal unitary operation U o (9 , 0i, m '_i) = I ] C k (9 k ). Now 

fc=0 

setting C/ O (0 O , 0i, m '-i) = G m {it) and p T (t) = b k and inserting them into 
the unitary transformation (38) one obtains 

d(m)— ld(m)— 1 

= b k -2[b k ,g m ]+ + A J2 J2 D k'bkD V) (39) 

k> v 

where D k >,Di> e S zq {m) x S zq (m). Since an anti-diagonal operator 6^ with 
an even index fc, e.g., fc m does not contain any diagonal operator component, 
that is, all the diagonal elements of the matrix b k are equal to zero, there must 
be D k ib k D k i = (b k ) k r k >D k i = for any diagonal operator D k i. On the other 
hand, there also holds: D k ib km Dy = (k' ^ I'), D k i,Dy e S zq (m) x S zq (m) 
(0 < m < n/2). The matrix element of D k /b k Di> is given by 

[D k ib k D v )ij = Sjk'dtk' {bk)ts&si'&ji' = 5ik'(bk)k'i'5ji>. 

t s 

The element is not zero only when i = k' and j = I' and the indices fc' and V 
satisfy the anti-diagonal line equation of the matrix b k :k' = —V + iV — 1 — k. 
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Since D k >, Dy e S zq {m) x S zq (m) the indices fc' and I' satisfy l m < k',l' < L m 
and k' + /' < 2L m . However, it follows from the anti-diagonal line equation 
that k' + 1' — N — 1 — k > 2L m because the condition (29) shows that the 
index N — 1 — k > 2L m . Therefore, the operator D k /b km Di/ is zero for those 
anti-diagonal operators b km used to build up the p— quantum operator Q pm . 
Then the unitary transformation (39) can be further reduced to the form 

GMhA^r 1 = b km - 2[b km ,g m ] + . (40) 

With the help of the unitary transformation (40) and the Trotter-Suzuki for- 
mula [13, 14] the subspace-selective p— quantum unitary operator U pm (9) = 
exp(—i0Q pm ) can be decomposed as 

U pm (6) = exp(-iehg m ,b km ]+) 

= exp[-i^6{b km - G m (7r)6 fem G' m (7r)" 1 )] 
= [exp(-i^6 fcm /L)G m (7r) 

x exp(i±8b km /L)G m (n)- 1 ] L + 0{L~ l ). (41) 

Note that norm for the operators b km and g m and their commutator [b km , g m ] 
equals one, that is, ||6 fe J| = 1, ||g- m || = 1, and || [b km , g m ] \\ = 1. Then for 
a modest number L = e^ 1 the expansion (41) for the p— quantum unitary 
operator U pm {6) can converge quickly. The computational complexity for 
the quantum circuit of the p— quantum unitary operator U pm (6) is therefore 
dependent on that of the unitary operations B km (6) and G m (9). It is shown 
in section 3.1 that the unitary operation G m (9) can be decomposed into a 
sequence of 2n basic unitary operations at most. For the situation that an 
unknown state in a subspace is transferred into the largest subspace of the 
Hilbert space of the n— qubit spin system the complexity of quantum circuit 
of the unitary operator B km (9) is ~ 0(log2 n), and for a general case that an 
unknown state is transferred from a subspace into a larger subspace the com- 
plexity is 0(2n 2 + Qne^ 1 ). Therefore, it follows from (41) that the subspace- 
selective p— quantum unitary operator U pm {9) can be expressed as a sequence 
of the basic unitary operations with complexity 0(2(2n 2 +6ne^ 1 )e^ 1 +Ans~ l ) . 

4. Discussion 
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It has been shown that any unknown quantum state can be efficiently 
transferred from a state subspace into a larger state subspace of the Hilbert 
space of an n— qubit spin system by a subspace-selective multiple-quantum 
unitary transformation, but the Grover quantum search algorithm [4, 15, 16] 
shows indirectly that the inverse process usually is hard one. This multiple- 
quantum transition process is similar to evolution process from nonequilib- 
rium state to equilibrium state in a closed quantum system, although the 
former is a unitary process and the latter a non-unitary and irreversible 
process [17]. This result might be helpful for understanding nonequilibrium 
processes such as protein folding process in nature where energy effect usually 
is not dominating from the point view of quantum dynamics. By the efficient 
subspace-selective multiple-quantum unitary transformation that can trans- 
fer efficiently any state from a small subspace into the largest subspace of the 
Hilbert space the search space of the quantum search problem can be reduced 
from the whole Hilbert space to its largest subspace. With the help of the re- 
sults in the paper and the auxiliary oracle unitary operation it can be shown 
that the quantum search algorithm [2, 6] based on quantum dynamics is at 
least as powerful as those quantum search algorithms including the Grover 
quantum search algorithm [4, 16] and adiabatic quantum search algorithm 
[18] in a pure quantum-state system because the former algorithm needs only 
to find which subspace the marked state is in among the (n + 1) subspaces 
of the Hilbert space. The diagonal and the anti-diagonal unitary operators 
will have an extensive application in constructing efficient quantum circuits 
for permutation operations of a symmetry group and the unitary operations 
of a cyclic group. 
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